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SUMMARY 


OBJECTIVE 

The  objectives  of  this  study  were  ( I )  to  determine  the  effect  of  various  REVGEN 
update  algorithms  on  the  output  of  a  matched  filter  processor,  and  (2)  to  evaluate  the 
resulting  false  alarm  and  detection  statistics  for  the  response  of  a  two-dimensional  variable 
threshold  detector,  which  operates  on  the  amplitude  range-Doppler  map. 

RESULTS 

REVGEN  algorithms  which  restrict  the  number  of  updates  and  randomize  in  phase 
only,  result  in  an  artificial  upper  limit  on  large  amplitude  peaks  in  the  range-Doppler  map. 
Thus,  they  yield  artificially  low  probabilities  of  false  alarm.  In  addition,  for  a  given  signal- 
to-noise  ratio,  these  lowered  rates  contribute  to  unreliably  high  probabilities  of  detec¬ 
tion.  due  to  incorrect  placement  of  detection  thresholds.  Overall,  it  is  seen  that  REVGEN 
algorithms,  which  generate  relatively  few  fixed-amplitude  updates,  can  have  a  significant 
effect  on  system  performance. 

RECOMMENDATION 

Both  theoretical  and  experimental  evidence  indicate  that  increasing  the  number  of 
updates  resolves  the  statistical  disparity.  More  significantly,  however,  randomizing  the  up¬ 
date  amplitudes  according  to  a  Rayleigh  distribution  produces  proper  single  bin  statistics 
independent  of  the  update  number.  This  decoupling  of  detection  statistics  from  update 
number  may,  in  some  instances,  permit  reduction  in  REVGEN  computational  requirements 
by  eliminating  the  need  to  process  a  large  number  of  updates.  Therefore,  the  randomiza¬ 
tion  algorithm  that  is  consistent  with  the  physics  of  the  medium  and  the  REVGEN  concept 
should  provide  tor  randomization  both  in  phase  and  in  amplitude. 
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1.0  INTRODUCTION 


In  the  rapidly  advancing  areas  of  undersea  weapon  technology  and  signal  processing, 
increased  emphasis  is  being  directed  toward  the  development  of  techniques  to  produce 
accurate  ocean  reverberation  signals.  The  use  of  computer  simulation  to  model  the  ocean 
acoustical  environment  can  reduce  both  the  time  and  expense  of  actual  at-sea  studies,  while 
providing  the  means  for  laboratory-based  investigations  essential  to  the  development  and 
testing  of  active  sonar  systems.  Much  effort  has  been  focused  on  real-time  implementation 
of  REVGEN  (Reverberation  Generator)  (ref.  1),  a  digital  simulation  of  ocean  reverberation. 

At  the  base  of  the  REVGEN  concept  is  the  assumption  that  reverberation  may  be 
represented  by  the  linear  process 

Y  =  B  X  (1) 

where  the  column  vector  X  is  a  sampled  version  of  the  transmit  pulse.  Y  is  a  sampled  version 
of  the  return  echo,  and  the  B  matrix  (a  “large”  set  of  complex  numbers)  contains  the  relevant 
information  about  the  ocean  environment  and  torpedo  ocean  interaction.  The  physical 
derivations  of  the  original  REVGEN  concept  are  documented  in  references  1 ,3,4,5  but.  as  a 
brief  outline,  can  be  summarized  in  the  following  three-step  process.  First,  a  Doppler 
density  matrix,  Dj  is  calculated  from  data  defining  the  ocean  environment  and  torpedo 
state  (i.e.,  velocity,  beam  patterns,  etc.)  which  in  turn  is  used  to  form  the  B-matrix  of 
equation  ( 1 ).  Finally,  the  B-matrix  and  samples  of  the  transmit  pulse  can  then  be  used  to 
generate  samples  of  a  simulated  reverberation  return. 

However,  the  requirement  that  this  ocean  reverberation  algorithm  operate  in  real 
time  on  the  present  NOSC,  Univac  1 1 10-based  Hybrid  Simulator  has  necessitated  several 
alterations  that  simplify  the  original  REVGEN  concept.  Although  an  accurate  REVGEN 
computation  would  typically  require  calculating  256  rows  of  256  Doppler  density  values 
every  pulse  length  (Tp),  the  current  capacity  of  the  simulator  is  to  calculate  one  range  slice  of 
256  Doppler  values  each  Tp.  This  disparity  between  the  full  number  of  range  slices  of  Dj  ^ 
and  the  fact  that  only  one  can  be  calculated  per  pulse  length  has  led  to  certain  “update” 
algorithms  which  serve  to  fill  parts  of  the  Doppler  density  matrix  without  resorting  to 
lengthy  calculation  of  the  ocean  algorithm.  The  ocean  algorithm,  used  to  construct  the  NOSC 
Doppler  density  matrix,  partitions  up  the  ocean  into  cells  according  to  range  and  Doppler 
and  then  calculates  the  deterministic  (mean)  energy  of  the  return  from  each  cell  according 
to  the  standard  sonar  equations  (ref.  4).  These  values  are  then  given  complex  random 


'  D.W.  Princehouse,  “REVGEN,  A  Real-Time  Reverberation  Generator  —  Concept  Development,”  Report 
No.  7511,  Applied  Physics  Laboratory,  University  of  Washington,  September  1975. 

“A.W.  Rihac/.ek.  “Principles  of  High  Resolution  Radar,”  McGraw-Hill,  New  York,  1969, 

~^D.W.  Princehouse.  "Reverberation  Generator  Ocean  Algorithm.  A  Status  Report,”  Report  No.  7806. 
Applied  Physics  Laboratory,  University  of  Washington.  February  1978. 

^B.A.  Bologna  and  E.M.  Rife,  “Real-Time  Ocean  Model  for  Reverberation,”  Naval  Ocean  Systems  Center, 
NOSC  TN  383,  March  1978. 

■’j.G.  Melville  and  M.E.  Stegman,  “The  Effect  of  Update  Randomization  on  REVGEN  Output:  Broadband 
Energy  Detector,”  Naval  Ocean  Systems  Center,  NOSC  TR  292,  August  1978. 


5 


phase  and  are  loaded  into  the  appropriate  element  of  the  matrix  Dj  j^.  This  departs  from 
the  true  form  of  the  matrix  in  two  ways:  first,  the  limited  update  algorithm  severely  reduces 
the  number  of  non-zero  matrix  elements;  second,  these  Dj  elements  are  assigned  deter¬ 
ministic  magnitudes  (calculated  from  the  sonar  equations),  when  in  reality  they  are  stochastic 
in  nature.  This  is  due  to  the  fact  that,  in  practice,  the  ocean  cells  contain  many  randomly 
placed  scatterers  whose  returns  add  coherently. 

This  report  examines  the  statistical  effect  of  these  modifications  (i.e.,  limited  updates, 
fixed  amplitudes)  on  detection  processing  of  the  simulated  signals.  (The  statistical  deviations 
addressed  here,  are  similar  to  those  presented  in  reference  5  for  the  response  of  a  broadband 
energy  detector.) 

By  analyzing  the  output  range-Doppler  map  of  a  matched  filter  processor*,  the  re¬ 
mainder  of  this  report  examines  the  effect  of  simplifications  of  the  REVGEN  algorithm  on 
corresponding  probabilities  of  false  alarm  and  detection.  Further,  where  changes  in  perform¬ 
ance  are  significant,  an  alternate  algorithm  is  presented  which  is  compatible  with  real  time 
operation  on  the  NOSC  Hybrid  Simulator  and  can  be  used  to  improve  system  performance. 


2.0  THEORY 

This  section  outlines  the  statistical  nature  of  (1)  the  output  obtained  from  matched 
filter  processing  of  reverberation  from  point  scatterers  in  an  ideal  ocean,  (2)  the  correspond¬ 
ing  analysis  applied  to  REVGEN  simulated  returns,  and  (3)  the  resulting  probabilities  of 
false  alarm  and  detection  for  the  response  of  a  two-dimensional  variable  threshold  detector 
which  operates  on  the  amplitude  range-Doppler  map. 


2. 1  RANGE-DOPPLER  MAP  FOR  POINT  SCATTERER  MODEL 

The  Generalized  Point  Scatterer  Model,  the  physical  model  from  which  REVGEN  was 
developed,  assumes  that  reverberation  is  the  sum  of  multiple  reflections  of  the  transmit 
pulse  from  point  scatterers  in  the  ocean  medium.  Each  echo  represents  an  individually 
amplitude-scaled,  time-delayed,  and  Doppler  shifted  version  of  the  transmitted  waveform 
(ref.  1 ).  This  can  be  expressed  mathematically  (ref.  3)  as: 


Y(t)  =  ^  A;  x  [tft-n)]  ;  (2) 

i=l  L  J 

where, 

Y  =  reverberation  signal  (baseband  complex) 

X  =  transmitted  signal  (baseband  complex) 

Aj  =  complex  amplitude  of  the  i1*1  scatterer,  including  scatterer  cross  section, 
beam  patterns,  attenuation,  propagation  loss,  and  phase  shift,  if  any, 
associated  with  reflection. 

‘Throughout  this  report,  “matched  filter”  processing  is  used  to  refer  to  the  operation  of  a  “replicate 
correlator,”  which  actually  “matches”  a  replica  of  a  transmit  pulse  to  the  return  signal. 
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7j  =  Doppler  shift  of  i1*1  scatterer  (approximately  1  +  2Vj/C’) 

Vj  =  speed  of  scatterer  relative  to  transducer  receiver 
C  =  speed  of  sound,  in  water 
Tj  =  (2  Rj/C)  =  time  delay  of  i1*1  scatterer 
Rj  =  range  of  i1*1  scatterer,  metres 
Ns  =  total  number  of  scatterers  in  ocean,  a  very  large  number. 

For  low  Doppler  targets  and  narrowband  transmissions  this  can  be  rewritten  in  the  approxi¬ 
mate  form 


Y(t)  =  ^  Aj  X(t-Tj)cjWi  <t_ti>  :  (3) 

i=  1 

where,  toj  is  the  received-frequency-caused  Doppler  shifts  of  the  individual  scatterers. 

Applying  these  results  to  the  formula  for  x(r,w),  (ref.  2)  the  matched  filter  (replicate 
correlator)  response  in  delay  and  Doppler  gives: 


OO 

XXy(r,co)  =  j  X(t)  Y*  (t-r)  el'wtdt 

— -OO 


(4) 


,  *  .  j(w-C0:)t  jcodr+tj) 

X(t)  x  (t-tj-r)  e  1  e  1  1  dt 


V  *  jwdr+tj)  ,r 

Xxy(T,co)  =  AiXxx<T  +  tj,w-0Ji)e 

i 

where  xxx  is  the  ambiguity  function  of  the  transmit  pulse. 

Thus,  the  output  of  the  matched  filter  processor  can  be  seen  to  equal  the  sum  of 
multiple  complex  scatterer  amplitudes  which  have  been  convolved  with  the  ambiguity  func¬ 
tion  of  the  transmit  pulse,  plus  a  phase  term. 

If  the  width  of  the  main  lobe  of  x(t,co)  is  greater  than  a  few  wave  lengths  of  the 
carrier  frequency  (which  is  almost  always  the  case)  and  when  there  arc  a  large  number  of 
scatterers  within  the  “main  lobe”  then  equation  (5)  represents  the  sum  of  a  large  number  of 
randomly  phased  components.  As  such,  x(t,co)  is  a  complex  Gaussian  random  variable  (by 
the  law  of  large  numbers)  and  lx(Y.co)|  has  a  Rayleigh  distribution.  Ideally  then,  for  the 
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point  scatterer  model  of  an  ocean  medium,  the  amplitudes  of  the  output  bins  in  the  range- 
Doppler  map  will  be  distributed  according  to  a  Rayleigh  distribution  (ref.  6). 

The  foregoing  has  reviewed  the  statistical  characteristics  of  the  range-Doppler  map 
produced  from  reverberation  that  would  be  received  in  an  idealized  (i.e.,  a  homogeneous 
medium  embedded  with  point  scatterers),  but  real  ocean  environment.  The  following 
section  discusses  the  effect  of  REVGEN  simulation  on  corresponding  matched  filter  output. 


2.2  RANGE-DOPPLER  MAP  FOR  REVGEN 

The  original  REVGEN  concept  has  been  thoroughly  documented  in  references  1,3 
but  is  briefly  outlined  here  to  establish  the  basis  of  the  relationship  between  the  elements  in 
the  computer  simulation  and  the  resulting  output  of  a  matched  filter  processor. 

Assuming  that  the  ocean  is  a  linear  medium  and  that  pulse  compression  or  expansion 
can  be  approximated  by  frequency  translation,  the  output  REVGEN  algorithm  can  be 
written  in  terms  of  the  following  matrix  equation: 

YJ=SBJXXL  (6) 

L 

where  represents  a  sampled  vector  of  the  transmit  pulse  and  Yj  correspondingly  represents 
the  return  signal.  If  and  Yj  are  samples  of  the  time  series  taken  at  intervals  equal  to  the 
inverse  of  twice  the  respective  bandwidths  (tox,coy)  of  the  transmit  pulse  and  return  echo,  it 
can  be  shown  that  the  B  matrix  contains  all  the  information  concerning  the  ocean  medium  that 
is  available  to  that  specific  ping-echo  process  (ref.  1 ). 

The  B  matrix  is  calculated  +  from  the  Doppler  density  matrix  D,  such  that 


N-l 

“J.L  "  £ 


M=0 


Dj-L.M<2’iML,N  = 


(7) 


where,  N  is  the  number  or  transmit  pulse  samples  in  X,(Lm^N-l),  and  J  indexes  the  time  of 
the  echo  return  (ref.  1 ).  The  Doppler  density  matrix,  Djjyj,  has  elements  that  correspond  to 
different  range  and  Doppler  returns  in  the  echo.  The  rows  of  Djj^  represent  the  spectral 
distribution  of  the  scatterers  at  a  particular  range  slice  J.  Throughout  this  report,  the 
Doppler  quantization  or  bin  size  is  taken  to  be  equal  to  1/Tp.  This  corresponds  to  the  finest 
Doppler  quantization  available  to  a  matched  filter  processor  for  a  ping  with  that  given  length. 

Given  the  Doppler  density  matrix,  REVGEN  then  produces  a  complex  time  series 
(demodulated)  which  corresponds  to  the  reverberation  return.  The  elements  of  the  Doppler 
density  matrix,  i.e.. 


DJ.K 


(8) 


^Rayleigh,  “Scientific  Paper,”  Vol.  VI, Cambridge  University  Press,  London,  1920. 

+The  sign  in  the  exponent  is  somewhat  arbitrary,  the  present  convention  is  consistent  with  positive 
frequencies  of  the  form  eltJt. 
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are  computed  (ref.  3)  from  the  range-  and  Doppler-quantized,  point  scatterers  of  the 
Generalized  Point  Scatterer  Model,  equation  (1),  such  that 


D 


J.K 


'J.K 


rJ0J,K 


I 


Aj  e 
certain  1 


-j0i 


(9) 


where,  the  sum  is  taken  over  those  scatterers,  i,  which  fall  into  the  range  cell  J  and  Doppler 
cell  K  and  where  4>-{  (and  therefore  the  cumulative  phase  0j  j^)  is  a  uniformly  distributed 
random  variable  between  zero  and  27r.  The  effect  of  this  quantization  is  discussed  in 
reference  1  and  is  shown  to  be  small. 

Applying  this  result  to  equation  (3)  for  the  output  of  the  matched  filter  processor, 
it  is  seen  that  for  REVGEN  returns  the  range-DoppIer  map  will  be  determined  by  the  sum  f 
the  amplitudes  in  the  Doppler  density  matrix  convolved  with  the  ambiguity  function  of  the 
transmit  pulse.  That  is. 


X(r.cj) 


X  °J,K  X<r+Tj  -  w-cok)eJ°Jk(r+TJ) 
J.K. 


(10) 


Using  equation  (9)  and  the  fact  that,  independent  of  the  distribution  of  any  variable  6. 
ej0+O  jlas  uniform  random  phase  if  eJ0  is  uniformly  distributed  in  phase,  it  can  be  shown 
that  equation  ( 10)  is  equivalent  to  the  randomly  phased  sum  of  the  products  of  the  magnitude 
of  Dj  and  the  magnitude  of  the  ambiguity  function  at  the  location  specified  by  the  J's  and 
K’s. 

The  following  section  compares  this  form  of  the  equation  with  the  original  theory 
given  in  equation  (5)  and  likewise  compares  the  probability  distributions  for  matched  filter 
processing  on  the  resulting  time  series. 


2.3  EFFECT  OF  REVGEN  MODIFICATIONS  ON  DETECTION  STATISTICS 

In  the  previous  sections,  two  expressions  were  developed  for  the  range-Doppler  maps 
from  an  idealized  point  scattering  model  and  for  a  range  and  Doppler  quantized  REVGEN 
calculation.  While  equations  (5)  and  (10)  are  similar  in  form  they  can  produce  significantly 
differing  returns.  This  section  will  explore  some  of  these  differences. 

The  primary  difference  between  equations  (5)  and  ( 10)  is  in  the  number  of  terms  in 
the  sum.  The  great  reduction  in  number  of  terms,  which  comes  with  the  REVGEN  calcula¬ 
tion,  has  its  basis  in  equation  (9)  which  shows  how  a  large  number  of  physical  scatterers  can 
be  mapped  into  a  single  range-Doppler  cell  Dj  When  this  sum  is  performed,  the  Aj 
values  (for  a  given  J  and  K)  should  be  distributed  according  to  a  Rayleigh  distribution,  with 
a  mean  value  which  corresponds  to  the  mean  return  from  the  portion  of  the  ocean  correspond¬ 
ing  to  the  J.K  range-Doppler  interval.  The  mean  returns  can  be  calculated  from  the  standard 
sonar  equations  and  a  knowledge  of  the  torpedo  beam  patterns. 

Up  to  this  point,  output  of  the  point  scattering  model  and  the  REVGEN  calculation 
should  be  identical  to  within  the  processing  constraints  imposed  by  the  input  and  output 
signal  bandwidths. 
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The  constraint  that  the  entire  reverberation  calculation  run  in  real  time  requires 
certain  reductions  to  the  original  REVGEN  algorithm.  In  particular,  the  number  of  times 
per  pulse  length  that  the  ocean  algorithm  can  be  executed  is  very  limited  (currently  only 
once  per  pulse  length).  Neither  the  point  based  method  of  equation  (9)  (which  gives  the 
proper  randomness  to  Dj  ^ )  nor  even  the  deterministic  calculation  of  Dj  from  the  sonar 
equations  can  be  calculated  at  the  required  rate  (which  can  be  up  to  hundreds  of  times  per 
pulse  length).  _  _____ 

The  current  approximation  make  use  of  the  fact  that  Dj  ^  *  Dj+r  K  (i.e.,  the  mean 

value  of  return  averaged  over  many  pings  changes  little  over  one  ping  length).  Each  ping 
length  Dj  ^  is  calculated  using  the  standard  sonar  equations  for  the  range  slice  J  and  for  all 
Doppler  bins  K.  To  fill  in  the  Doppler  Density  matrix  at  update  locations  within  one 
ping  length,  rp  of  J  the  following  approximation  is  made  Dj+Aj  k  =  Dj  K' 

Randomness  is  introduced  for  each  new  update  by  multiplying  terms  Djj^  by 
gj^J+AJ.K  where  0J+AJ,K  is  selected  randomly  from  a  uniform  distribution.  The  values 
Dj+AJ  K  are  set  e9ua'  to  zero  f°r  AJ  =  1  to  (N§^p/N)-1 ,  where  N<jAmP  is  the  number  of 

range  samples  per  ping  length  and  N  is  the  number  of  updates  per  pulse  length..  The  range 
slice  J+NgAMp/N *s  then  created,  according  to  the  relationship 


‘>L.K  =  Dj.K«i*L,K 

where 

L  =  J+NSAMP/N- 


This  procedure  of  re-randomizing  (in  phase  only)  the  initial  range  slice  and  using  it  to 
“update”  the  Doppler-density  matrix  continues  at  regular  intervals  (^SAMp/hO  throughout 
the  ping. 

This  undersampling  of  the  Doppler-density  further  reduces  the  number  of  terms  in 
equation  (10)  and  can  directly  affect  the  statistics  for  the  magnitude  lx(r,to)|,  of  the  output 
of  a  matched  filter  processor.  When  the  Dj  K  are  fixed  according  to  a  deterministic  calcula¬ 
tion  of  the  standard  sonar  equations  (ref.  4,  based  on  the  mean  intensity  of  scatterers  from 
within  a  certain  volume)  the  sum  x(r,co)  varies  only  according  to  phase  angles  and  the  shading 
applied  by  the  ambiguity  function.  Therefore,  equation  (10)  has  an  absolute  limit  for  the 

possible  values  of  lx(r,cc)lmax  =  £  I  Dj-k*  I . 

Considering  the  probability  distribution,  it  is  consequently  impossible  to  obtain 
arbitrarily  large  values  of  lx(r,w)|  beyond  this  implied  maximum.  (This  is  identical  to  the 
problem  discussed  in  reference  5  for  the  envelope  statistics  of  the  complex  time  series 
generated  by  REVGEN.) 

When  the  number  of  significant  components  in  the  sum  (i.e.,  updates  times  the 
number  of  filled  bins  per  update  row)  is  large,  this  maximum  value  is  large  compared  to  the 
mean  value  and  the  distribution  is  more  realistic  in  the  tails;  that  is,  as  the  number  of  updates 
increases,  the  distribution  of  the  amplitudes  in  the  range-Doppler  map  will  approach  a 
Rayleigh  distribution.  However,  when  the  number  of  components  is  small  (<20),  significant 
differences  can  be  noted  in  the  shape  of  the  tail  of  the  distribution  -  the  area  crucial  in 
determining  false  alarm  rates.  Since  it  is  impossible  to  obtain  values  greater  than  the  previously 
defined  maximum,  the  probability  distribution  beyond  this  point  drops  sharply  to  zero. 
Applying  this  to  the  probability  of  false  alarm,  PpA>  which  is  determined  by  the  area  in  the 
tail  of  the  distribution  beyond  a  set  detection  threshold,  it  is  seen  that  as  the  placement  of 
the  detection  threshold  approaches  the  “drop-off’  point,  the  values  obtained  for  PpA  will 
be  erroneously  low. 
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However,  there  is  an  alternative  to  the  time  consuming  process  of  increasing 
the  number  of  updates  in  order  to  improve  the  statistics  of  REVGEN  simulation.  If  the 
elements  Dj  ^  of  the  Doppler  density  matrix  are  randomized  such  that  Aj  is  Rayleigh 
distributed  (with  the  mean  of  the  distribution  set  by  the  deterministic  sonar  equation 
calculations),  then  the  magnitudes  Ix(t,cj  )l  will  have  the  proper  statistical  distribution 
independent  of  the  number  of  elements  used  in  the  computation.  In  this  case,  equation  ( 10) 
will  be  the  sum  of  complex  Gaussian  random  variables  and  is,  therefore,  itself  a  complex 
Gaussian  random  variable  having  the  property  that  its  magnitude  Ix(r.cc)  I  follows  a 
Rayleigh  distribution. 

Therefore,  a  Rayleigh  distribution  of  the  magnitudes  of  the  matched  filter  output 
can  be  achieved  in  either  of  two  ways:  (1 )  by  using  many  updates,  or  (2)  by  randomizing 
the  magnitudes  of  the  Doppler  density  matrix  elements  according  to  a  Rayleigh  distribution. 
Thus,  the  REVGEN  update  algorithm  used  to  generate  the  elements  of  the  Doppler  density 
matrix  has  a  direct  effect  on  the  resulting  output  from  a  matched  filter  processor  and  can 
therefore  alter  the  response  of  a  two-dimensional  variable  threshold  detector  which  operates 
on  the  amplitude  range-Doppler  map. 


3.0  RESULTS 

Based  on  the  theories  just  presented,  the  results  of  this  section  provide  graphic- 
illustrations  and  experimental  evidence  of  the  effect  various  REVGEN  update  algorithms 
have  on  the  probability  statistics  for  the  output  from  a  matched  filter  processor.  By 
analyzing  the  distribution  of  the  amplitudes  in  the  output  bins  (separated  by  1  /r„  in 
Doppler)  of  the  range-Doppler  map,  the  following  figures  provide  a  basis  for  comparing  the 
various  false  alarm  and  detection  statistics  for  the  response  of  a  two-dimensional  variable 
threshold  detector. 

In  addition  to  the  predicted  theoretical  values,  results  are  also  demonstrated  with 
actual  REVGEN  data.  Computer  simulation  using  a  rectangular.  CW  transmit  pulse.  20 
channels  and  a  three-second  listening  period  generated  20  000  independent  samples  of  the 
output  range-Doppler  map.  To  reduce  the  time  and  cost  of  computation  it  was  determined 
that  every  eighth  bin  in  each  Doppler  density  row  could  be  filled  and  still  maintain  independ¬ 
ence  between  samples.  Generally,  it  would  be  expected  that  some  leakage  would  occur  in 
adjacent  bins.  However,  it  can  be  shown  that  during  simulations  with  2.  4.  and  8  updates 
per  pulse  length  when  the  Doppler  values  are  spaced  with  eight  bin  separation,  all  other 
non-zero  elements  of  a  row  will  not  contribute  the  resulting  matched  filter  output  because 
the  value  of  the  ambiguity  function  will,  in  fact,  be  zero  at  these  intervals.  (See  Appendix.) 


3.1  DISTRIBUTION  ANALYSIS 

In  the  final  theory  section,  emphasis  was  placed  on  the  comparison  of  Rayleigh 
statistics  for  an  idealized  ocean  model  and  the  corresponding  statistics  for  REVGEN  return 
generated  with  relatively  few,  fixed  amplitude,  uniform  phase  distributed  updates.  It  was 
shown  that  the  output  of  a  matched  filter  processor  is  significantly  affected  when  REVGEN's 
updates  are  limited  in  number  and  fixed  in  amplitude.  The  comparison  of  figures  1  and  2 
illustrates  the  effect  of  update  randomization  on  the  amplitude  distribution  of  the  matched 
filter  output  of  actual  REVGEN  data  which  was  generated  with  two  updates  per  pulse 
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length.  Figure  1  shows  the  distribution  for  updates  of  constant  amplitude  and  uniform 
random  phase.  Figure  2  is  the  experimental  result  for  a  filtered  REVGEN  return  when  the 
two  updates  per  pulse  length  were  given  Rayleigh  randomized  amplitudes.  (In  both  figures, 
the  solid  line  represents  theoretical  Rayleigh  statistics  with  the  same  sample  mean.  Error 
bars  represent  one  standard  deviation.) 


3.2  PROBABILITY  OF  FALSE  ALARM 

In  processing  simulated  reverberation  returns,  a  two-dimensional  variable  threshold 
detector  operates  on  the  amplitudes  in  the  output  bins  of  the  range-Doppler  map  and 
responds  with  a  “target”  detection  when  this  amplitude  exceeds  both  a  set  detection 
threshold  and  a  variable  threshold  based  on  the  average  amplitudes  in  surrounding  bins. 

The  altered  amplitude  statistics  for  Dj^  shown  in  the  previous  section  directly  affect  the 
probability  of  false  alarm  since  it  is  precisely  the  probability  of  obtaining  high  amplitude 
range-Doppler  peaks  (i.e.,  the  area  in  the  upper  tail  of  the  distribution)  which  is  critical  in 
determining  false  alarm  rates.  The  artificial  upper  limit  caused  by  REVGEN  algorithms 
which  use  relatively  few  updates,  not-randomized  in  amplitude,  results  in  artificially  low 
values  for  Pfa 

In  terms  of  equation  (10),  lxXy  (t,co)|  will  approach  a  Rayleigh  distribution  even 
for  updates  whose  magnitudes  are  constant  if  the  number  of  Dj  which  contribute 
significantly  to  the  sum,  is  large.  The  exact  weighting  of  each  component  is  determined  by 
ambiguity  function  xXx  (T  +  TJ>  00  +  UK)  For  a  square  pulse,  this  function  falls  off  rapidly 
in  Doppler  and  relatively  slowly  in  range  (unless  explicitly  noted  all  analysis  in  this  section 
refers  to  a  square  pulse  envelope,  the  effect  of  pulse  type  is  treated  in  the  discussion  section). 

Figures  3  and  4  compare  the  effect  on  the  false  alarm  rate  (for  a  fixed  threshold)  of 
increasing  the  number  of  significant  contributors  to  equation  (10).  This  can  be  done  in 
two  ways,  increasing  the  number  of  Doppler  bins  filled  per  range  slice  (figure  3)  and  by 
increasing  the  number  of  updates  per  pulse  length  (figure  4).  In  both  cases,  the  curves  were 
determined  by  first  defining  a  set  of  amplitudes  according  of  the  ambiguity  function  and 
the  positions  of  the  filled  bins  in  the  proposed  Doppler-density  matrices.  The  problem  of 
finding  the  distribution  of  the  resultant  of  the  vector  sum  of  a  finite  set  of  randomly 
phased  components  (each  having  a  specified  length)  was  then  solved  as  in  reference  5.  This 
technique  was  used  to  produce  all  the  theoretical  curves  in  the  present  report. 

Figure  3  shows  the  probability  of  false  alarm  curves  for  various  numbers  of  Doppler 
bins  filled  (centered  on  the  maximum  Doppler  response  of  the  ambiguity  function)  for 
eight  updates  per  pulse  length.  As  the  number  of  filled  bins  increases,  the  probability  of 
false  alarm  becomes  more  like  the  Rayleigh  curve.  However,  since  the  ambiguity  function 
ceases  with  increasing  distance  from  the  center  bin,  additional  terms  are  of  decreasing 
importance. 

Figure  4  shows  the  probability  of  false  alarm  for  various  numbers  of  updates  per 
pulse  length  when  one  bin  is  filled.  The  limits  to  the  response,  resulting  from  the  small 
number  components,  is  clearly  seen.  Note  that  even  for  eight  updates  per  pulse  length  and 
a  PpA  the  false  alarm  rate  is  more  than  an  order  of  magnitude  below  the  theoretical 

value. 

For  experimental  support  of  the  theories  just  presented,  figures  5  through  7  present 
distributions  of  lxxy  ( r.w  )  1  derived  from  20  000  samples  of  the  range-Doppler  maps 
produced  from  actual  REVGEN  returns  generated  with  2,  4,  and  8  updates  per  pulse 
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Figure  6.  Same  as  figure  5  with  4  updates  per  pulse  length,  constant 
amplitudes. 


PFA 


<R-2!/2(SIGMA-2) 


Figure  7.  Same  as  figure  5  with  8  updates  per  pulse,  constant  amplitudes. 


length  (respectively),  constant  amplitude,  a  rectangular  transmit  pulse,  uniform  phase 
distribution  and  one  filled  Doppler  bin.  These  computer  simulations  yielded  results  which 
are  consistent  with  the  theoretical  values  of  figure  4  ( 1  bin  filled)  and  also  substantiate  the 
predicted  improvement  caused  by  an  increase  in  the  update  algorithm.  (Error  bars 
represent  one  standard  deviation.) 

For  comparison,  figure  8  presents  a  distribution  derived  from  20  000  samples  of  a 
REVGEN  filtered  return  which  were  generated  like  those  of  figures  5  through  7,  except  that 
the  amplitudes  of  the  two  updates  per  pulse  length  were  randomized  according  to  a  Rayleigh 
distribution.  Clearly,  it  is  seen  that  while  constant  amplitudes  require  a  large  update 
algorithm  to  produce  the  proper  statistics,  Rayleigh  randomized  updates  follow  the  Rayleigh 
statistical  model  independent  of  the  number  of  updates  used  to  generate  the  signal. 

Table  1  provides  further  evidence  that  REVGEN  algorithms,  which  randomize  the 
elements  of  the  Doppler  density  matrix  in  phase  only,  can  cause  artificial  lowering  of 
observed  false  alarm  rates.  Up  to  this  point,  it  has  been  shown  both  theoretically  and 
experimentally  that  the  probability  of  false  alarms  is  significantly  effected  by  changes  in  the 
update  algorithm.  Table  1  summarizes  actual  test  results  of  false  detections  which  occurred 
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Figure  8.  Same  as  figure  5  except  the  update  algorithm  uses  2  updates  per  pulse 
length  having  a  Rayleigh  distribution. 


Table  1.  Post  detection  editing  false  alarm  rates. 


REVGEN 

False 

Update  Algorithm 

Alarms 

2  updates/pulse  length  Rayleigh 
randomized  amplitudes,  uniform  phase 

24  ±5 

8  updates/pulse  length  Rayleigh 
randomized  amplitudes,  uniform  phase 

26  ±  5 

4  updates/pulse  length  fixed  amplitudes, 
uniform  phase 

5  ±  3 

8  updates/pulse  length  fixed  amplitudes, 
uniform  phase 

9  ±  3 
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when  a  representative  two-dimensional  variable  threshold  detector  was  applied  to  the 
matched  filter  output  (i.e.,  range-Doppler  maps)  of  REVGEN  signals  generated  with  various 
update  algorithms.  In  this  case,  all  frequency  bins  were  filled  and  the  number  of  indepen¬ 
dent  beams  was  increased  until  a  significant  number  of  false  alarms  had  been  detected  for 
each  of  the  update  randomization  schemes. 

3.3  SNR  POWER  CURVES 

The  power  curves  shown  in  figures  9  through  1 1  are  used  to  compare  the  relative 
accuracy  of  values  obtained  for  the  probability  of  detection  (Pq)  for  differing  numbers  of 
updates  per  pulse  length.  Figures  9  through  1 1  set  the  detection  threshold  at  a  predetermined 
PpA  level  and  examine  the  resulting  Pq  values  for  various  signal-to-noise  ratios.  Given  the 
altered  Pp^  statistics  just  illustrated,  it  is  clear  that  REVGEN  update  algorithms  will  like¬ 
wise  alter  the  SNR  results,  since  these  inaccurately  low  Pp^  rates  are  those  used  to  determine 
the  placement  of  threshold.  Thus,  the  focus  of  figures  9  through  1 1  indicates  that  at  various 
SNR  levels,  the  erroneously  low  PpA  values  produce  erroneously  high  P[>  rates  due  to  in¬ 
correct  placement  of  the  threshold. 

The  severity  of  this  problem  is  somewhat  dependent  upon  how  the  data  is  used.  If 
the  detection  threshold  is  determined  by  taking  a  local  average  of  values  in  the  range-Doppler 
map  and  then  setting  the  threshold  based  on  the  assumption  of  a  true  Rayleigh  distribution, 
then  the  threshold  will  be  close  to  the  correct  value  and  accurate  values  of  Pj)  will  be  meas¬ 
ured.  However,  in  this  case  the  PpA  values  will  be  too  low.  If  the  data  is  used  to  experi¬ 
mentally  set  the  detection  threshold  by  measuring  the  false  alarm  rate,  then  the  resulting 
Pq  will  be  too  high. 


19 


SNR  PLOT  SNR  PLOT 


Figure  9.  Probability  of  detection  as  a  function  of  SNR  Figure  10.  Same  as  figure  9  with  PFA  =  10~: 

(threshold  set  at  PFA  =  10"^)  for  various  update 
algorithms  and  for  the  Rayleigh  case. 


SNR  PLOT 


SNR  fdB) 

Figure  1 1.  Same  as  figure  9  with  PFA  =  10“^. 


4.0  DISCUSSION 

From  the  preceding  sections,  it  is  seen  that  the  update  algorithm  used  to  generate 
REVGEN  reverberation  returns  can  directly  affect  the  results  of  detection  processing  of  the 
simulated  signals.  Both  theoretical  and  experimental  analysis  of  the  output  of  a  matched 
filter  processor  showed  notably  altered  statistics  for  the  probabilities  of  false  alarm  and  de¬ 
tection  when  changes  were  made  in  the  update  algorithm. 

Most  significantly,  it  was  shown  that  simulations  which  limit  the  number  of  updates 
per  pulse  length,  while  randomizing  in  phase  only,  yield  artificially  low  Pp^  rates  and.  there¬ 
fore,  affect  proper  placement  of  detection  thresholds. 

While  increasing  the  number  of  updates  seems  to  resolve  the  disparity,  it  was  also 
shown  that  randomizing  the  magnitudes  according  to  a  Rayleigh  distribution  produces  the 
proper  distribution  of  values  within  a  single  range-Doppler  bin,  independent  of  the  number  of 
updates  used  in  the  simulation. 


This  decoupling  of  the  detection  statistics  fro...  the  update  algorithm  may,  in  some 
instances,  permit  a  reduction  in  the  computational  time  needed  for  real-time  implementation 
of  REVGEN,  by  eliminating  the  need  to  process  a  large  number  of  updates. 

While,  by  using  the  proper  randomization,  the  single  bin  values  can  be  made  to  have 
the  correct  distribution  with  few  updates  per  pulse  length,  there  are  other  features  of  the 
range-Doppler  map  which  remain  sensitive  to  the  number  of  updates  used.  These  effects  deal 
with  bin-to-bin  correlation  in  range  and  the  relative  mean  value  of  the  bins  from  one  range  to 
the  next.  A  brief  analysis  is  given  to  highlight  these  effects.  For  the  numerical  purposes  of 
this  analysis  a  rectangular  transmit  pulse  envelope  is  assumed.  This  is  noted  where  other 
pulse  types  would  make  a  qualitative  difference  in  the  results.  First,  the  bin-to-bin  correla¬ 
tion  is  analyzed. 

In  mathematical  terms  the  range-Doppler  map  results  from  the  two-dimensional 
convolution  of  the  ambiguity  function  of  the  transmit  pulse  with  the  Doppler  density 
matrix,  which  for  a  point  scattering  model,  is  a  Finite  set  of  delta  functions  distributed  in 
range  and  Doppler  (as  in  equation  (5)).  This  ambiguity  function  has  a  value  of  one  at  its 
center  and  a  non-zero  value  for  a  range  extent  of  two  pulse  lengths. 

As  seen  in  equation  (5),  the  sum  x(r\cc)  will  include  some  of  the  same  Aj’s  as  the  sum 
X(r+Ar,cu)  (if  |Arl  <  rp).  Therefore,  one  would  expect  x(T><^)  to  be  correlated  with 
X(r+Ar. co)  for  |Ar|  <  rp.  In  terms  of  a  matched  Filter  processor  of  the  GCCD-type  a  bin- 
to-bin  correlation  in  range  is  expected  for  range  slices  less  than  one  pulse  length  apart. 

It  does  not  follow,  that  the  correlation  for  a  given  overlap  (1  -Ar/rp)  is  independent 
of  the  number  of  updates.  A  few  computer  experiments  were  performed  to  demonstrate  this. 

For  this  case,  of  one  Doppler  bin  Filled  with  complex  Gaussian  random  numbers  (Aj), 
the  value  in  the  range-Doppler  map  xj  at  the  given  Doppler  has  the  form 


j=-(N-l ) 


Aj+j  R(J) 


where,  N  is  the  number  of  updates  per  pulse  length,  R(J)  is  the  autocorrelation  function  of 
the  transmit  pulse  and  it  is  assumed  that  the  xi  is  measured  at  the  range  corresponding  to 
Aj.  Equation  (12)  uses  a  single  subscript  for  x  to  denote  the  fact  that  only  values  at  one 
particular  Doppler  are  being  calculated  (the  Doppler  corresponding  to  the  one  filled  bin). 

For  a  rectangular  transmit  pulse  R(J)  =  |N-J|/N,  where  J  =  -N  to  N  and  otherwise 
zero;  therefore. 


N-l 

T  I 

j=-N-l 


vi+j  T 


Several  series  of  Xj  were  generated  from  random  sequences  of  Aj’s  using  equation  (13) 
with  successive  values  of  the  series  corresponding  to  a  change  in  the  observation  point  of  half 
the  pulse  length  (i.e.,  Ai  =  N/2).  This  was  done  for  N  =  8,  and  N  =  2. 
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The  autocorrelation  R(t)  function  was  then  calculated  for  the  function  Ixjl  of  the 
various  series  at  the  point  r  =  rp/2.  The  results  of  several  series  having  the  same  N,  but 
different  random  number  sequences  (Aj),  were  averaged  and  appear  in  table  2;  also,  included 
in  this  table  are  the  corresponding  correlations  taken  from  the  same  REVGEN  simulations 
used  to  produce  the  distributions  reported  earlier.  Also,  included  for  comparison  are  the 
correlation  values  for  the  REVGEN  runs  which  randomized  in  phase  only.  The  statistics  of 
this  problem  are  slightly  different  due  to  the  randomization  technique  and  so  no  predicted 
value  is  given. 

A  final  entry  in  table  2  is  the  expected  correlation  when  there  is  a  very  large  number 
of  updates.  The  calculation  is  similar  to  that  described  in  chapter  4  of  OTShevskii  (refer¬ 
ence  7),  except  that  the  “signal  envelope”  in  this  context  must  refer  to  the  ambiguity  func¬ 
tion  rather  than  the  envelope  of  the  transmit  pulse,  as  is  the  case  when  the  correlation  of  the 
reverberation  envelope  function  is  sought.  It  can  be  shown  that  both  problems  have  the 
same  mathematical  form,  if  the  above  substitution  is  made. 

Table  2.  Expected  and  observed  correlation  for  matched 
filter  outputs  which  are  separated  in  range  by  half  (square 
transmit  pulse). 


Number  of 
Updates/Pulse 

Magnitude 

Randomization 

R(Tp/2) 

from 

REVGEN 

R(rp/2) 

from 

Eq. (12) 

Theoretical 

(Ol’Shevskii) 

2 

Rayleigh 

.435  +  .007 

.416  +  .009 

8 

Rayleigh 

.495  +  .05 

.487  +  .021 

Constant 

,469  +  .007 

8 

Constant 

.492  +  .007 

Large 

Rayleigh 

.488 

There  are  two  general  observations  to  be  made  from  the  data  in  table  2.  The  first,  is 
that  as  the  number  of  updates  per  pulse  length  is  decreased  the  correlation  is  decreased 
between  two  bins  in  the  range-Doppler  map,  separated  in  range  by  half  a  pulse  length.  The 
second,  is  that  the  theoretical  value  for  large  N  is  within  the  expected  uncertainty  of  the 
correlation  for  eight  updates  per  pulse  length. 

For  shaped  transmit  pulses  the  ambiguity  function  falls  off  more  rapidly  in  range  than 
for  a  square  pulse.  This  results  in  small  weighting  for  those  scatterers  which  would  be  common 
to  two  range-Doppler  map  values,  separated  by  half  a  pulse  length  and  consequently  in  a  low 
bin-to-bin  correlation  for  even  large  values  of  N. 

The  other  general  feature  of  the  range-Doppler  map,  which  changes  with  update 
number,  is  the  expected  mean  value  of  Ixjl  at  ranges  between  those  range  slices  where 
scatterers  were  located  in  the  Doppler  density  matrix.  This  effect  can  most  easily  be  exam¬ 
ined  by  again  considering  the  case  of  a  Doppler  density  matrix  where  only  one  bin  is  filled 
in  Doppler  and  a  square  transmit  pulse  is  used. 


V  V.  Ol’Shevskii,  “Characteristics  of  Sea  Reverberation,”  Translated  from  Russian  by  Consultants  Bureau, 
New  York  1967. 
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If  the  scatterers  have  a  complex  Gaussian  distribution  with  a  =  Oj  then  the  expected 
value  of  lx(r)|  is  given  by 


lx(r)l 


^  (a j  R(j+T))2 
j=N-l 


(14) 


For  the  case  N=2,  oj=  1 ,  and  a  rectangular  transmit  pulse  with  R(J)  is  a  given  in  equation  (13). 


IX(0)I  = 


+  (la)2  + 


Vi 


=  1.22a 


Here,  lx(0)|  indicates  that  1x1  is  estimated  at  a  point  in  the  range-Doppler  where  one  of  the 
scatterers  existed.  In  similar  fashion,  bdTp/N)|  is  the  estimated  value  of  Ixl  one  update  length 
away  from  the  reference  position  (for  a  stationary  process,  with  all  propagation  losses  re¬ 
moved,  one  would  expect  lx(0)l  =  lx(mrp/N)|  where  m  is  an  integer.) 

In  the  above  example  the  weightings  of  1/2,  1,  and  1/2  result  from  the  form  of  the 
ambiguity  function  for  a  rectangular  pulse.  This  is  illustrated  in  figure  12. 

If,  however,  Ixl  is  to  be  estimated  at  a  point  half  way  in  between  scatterer  locations, 
the  shad.ng  of  the  ambiguity  function  (for  N=2)  is  1/4,  3/4,  3/4,  and  1/4;  therefore, 


lx(Tp/4)| 


=  [(H~+(t°)  Wt°)  +(?°)~  =  112°- 


At  Tpl 2,  the  shading  is  the  same  as  for  x(0)  and  therefore  lx(0)|  =  lx(rp/2)|. 

As  the  number  of  updates  per  pulse  length  increases,  the  difference  decreases  between 
the  maximum  value  (where  the  observation  is  at  the  same  range  slice  as  one  of  the  scatters) 
and  the  minimum  expected  value  (which  occurs  at  a  point  half  way  between  the  scatters). 

A  few  sample  curves  showing  this  variation  are  shown  in  figure  13. 

This  figure  shows  curves  for  various  pulse  types  and  update  numbers.  For  cosine 
and  cosine  squared  weighted  pulses,  the  following  formulas  were  used  for  R(r)  (from 
reference  7): 

for 

S(t)  =  cos(7rt/rp) ,  |t|<rp/2 


R(t)  =  g^-Jy-JcosOrr/rp)  + -^  sin  (7r|r|/Tp)j  (15) 

^ef  ~  rp/2 
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O 

b) 

0 

0 

c)  X(o) 

0 

0 

d)  X(Tp/4) 

0 


-2 


-1 


Figure  1 2.  a)  Aw  =  0  slice  of  ambiguity  function  for  square 
pulse  of  length  fp  b)  one  dimensional  range  doppler  map 
c)  representation  of  weighting  function  for  x(0)  d)  representa¬ 
tion  of  weighting  function  for  xOp/4)- 


25 


Figure  13.  Curves  of  lx(r)l/IX(0)|  for  various  update  rates 
and  pulse  types.  S2  =  square  pulse,  2  updates,  S4  =  square 
pulse  4  updates,  S8  =  square  pulse,  8  updates,  C2  =  cosine 
weighted  pulse,  2  updates,  CS2  =  cosine  squared  weighted 
pulse,  2  updates. 


and  for 


S(t)  =  cos20rt/Tn),  |t |  <  tJ2 


t^sln(2»Wfrp)-^sin(^l!)J 


(16; 


5  -3tP 

O-f - 


ef  8 


In  general,  the  modulation  in  the  amplitude  of  Ixl  is  slight.  The  worst  case  shown  in 
the  Figure  (two  updates  per  pulse  length  with  as  cos2  transmit  pulse)  has  a  maximum  devia¬ 
tion  of  x2  of  0.8  dB. 


Table  3  contains  numbers  representing  Ix(t0I/Ix(0)I  for  several  updated  rates  and 
pulse  types  where  r  was  chosen  so  that  Ix(t)|  is  at  a  minimum.  This  again  shows  that,  except 
for  the  case  of  one  update  per  pulse  length  (shown  separately  in  figure  14),  the  variation  in 
1x(t)1  is  small. 

Table  3.  Values  of  lx(rmin)l/lx(0)l  for  various  updates 
per  pulse  length  and  pulse  types.  S  =  square  pulse, 
c  =  cosine  weighted  pulse,  cs  =  cosine  squared  pulse. 


Updates 

Pulse 

Tmin 

1 

S 

V2 

0.707 

1 

C 

V2 

0.450 

1 

CS 

V2 

0.236 

2 

s 

V4 

0.913 

2 

c 

Tp/4 

0.976 

2 

cs 

Tp/4 

0.907 

4 

s 

t/8 

0.977 

4 

c 

t/8 

0.999 

4 

cs 

t/8 

0.999 

8 

s 

t/8 

0.994 

8 

c 

t/8 

1.000 

8 

cs 

t/8 

1.000 

The  severity  of  this  effect  on  a  simulation  will  depend,  in  part,  on  how  the  range- 
Doppler  map  is  sampled.  If  the  difference  in  range  between  samples  of  the  range-Doppler 
map  (the  FFT  overlap  is  a  matched  filter  type  processor)  is  an  integral  multiple  of  the 
distance  between  updates  in  the  Doppler  density  matrix,  then  the  bin  statistics  should  be 
stationary  throughout  a  ping.  Even  in  this  case,  small  variations  in  the  mean  level  of 
reverberation  can  occur  from  ping-to-ping,  depending  on  where  the  initial  sampling  interval 
is  placed  with  respect  to  the  update  scatterers.  When  the  sampling  frequency  is  not  an 
integral  divisor  of  the  update  rate  the  local  mean  reverberation  level  will  vary  periodically 
during  the  ping  cycle. 

Therefore,  it  is  seen  that  one  cannot  arbitrarily  reduce  the  update  rate  without  affect¬ 
ing  the  nature  of  the  REVGEN  output.  While  the  distribution  of  values  in  a  single  bin  of  the 
range-Doppler  map  is  unaffected  by  update  rate  or  pulse  type,  if  a  Rayleigh  randomization  is 
applied  to  the  updates,  the  bin  to  bin  correlation  in  range  and  the  mean  power  per  bin  as  a 
function  o Grange  are  still  affected.  These  last  two  factors  then  (upon  applying  the  proper 
update  randomization),  become  the  critical  issues  in  deciding  upon  the  proper  number  of 
updates  to  use  for  simulation  studies  with  matched  filter  type  processors. 


T/Tp 


Figure  14.  Curves  of  I X(r)  1/1  x(0)l  for  one  update  per  pulse 
length  and  various  pulse  types.  S  =  square  pulse,  C  =  cosine 
weighted  pulse,  CS  =  cosine  squared  weighted  pulse. 
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5.0  CONCLUSION 


The  conclusion  to  be  drawn  from  the  preceding  paragraphs  is  that  the  REVGEN 
algorithm,  which  determines  the  number  and  randomization  of  the  updates,  has  a  direct 
effect  on  the  output  of  a  matched  filter  processor.  This,  in  turn,  can  significantly  affect  the 
false  alarm  and  detection  statistics  for  target  detectors  which  operate  on  the  amplitude 
range-Doppler  map. 

Specifically,  when  the  updates  are  limited  (due  to  real-time  restrictions)  and  random¬ 
ized  in  phase  only,  the  resulting  Pp^  statistics  are  artificially  low;  thus,  forcing  the  placement 
of  detection  thresholds  to  be  correspondingly  inaccurate.  Both  theoretical  and  simulated 
evidence  indicate  that  increasing  the  number  of  updates  resolves  the  disparity.  However,  it 
is  significant  to  note  that  Rayleigh  randomized  updates  produce  the  proper  single  bin 
statistics,  independent  of  the  number  of  updates  used  in  the  simulation.  Certain  bin-to-bin 
characteristics  are  affected  by  reduction  in  update  number,  but  these  are  small  and  the 
effects  are  well  understood.  The  decoupling  of  the  detection  statistics  from  the  update 
algorithm  may,  in  some  instances,  permit  reduction  in  REVGEN  computational  requirements 
by  obviating  the  need  to  process  a  large  number  of  updates  per  pulse  length. 
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APPENDIX 


To  determine  the  form  of  the  ambiguity  function,  consider  the  baseband  equation 

(ref.  2) 

OO 

X(t,o>)  =  J  x(t)  x*  (t  -  r)  eiu,t  dt  ; 

— OO 


io>  t 

where  x(t)  is  of  the  form  x(t)  =  e  °  rect  (t). 

Substitution  yields: 

OO 

X(t,cj)  =  J  e  W°  rect  (t)  e  rect(t -t)  dt 

— OO 


OO 

=  J  e*^0*  rect  (t)  e  rect(t-T)  e‘wt  dt 

— OO 


OO 

=  elW°T  J  rect(t)  rect(t-r)  e’wt  dt  ; 

— OO 


which,  integrated  over  the  interval  T'  =  T  -  r  gives 


Writing  cj  in  the  form  cj  =  (n/T)2jr,  where  wn  is  the  frequency  in  the  nth  Doppler  bin,  we 
can  substitute 


X(r,w)  = 


_  sin  ai(T-r)/2 


cu/2 


X(r,n)  = 


sin  y  2rr(T -r)/ 2  sin  nff  ^  1  -  y) 


nff 

T 


(T-t)  sin  nff  ^1 


(T-r) 


nrr 


=  (T-r) 


sin  nn  j 

t 

t) 

ntr  ) 

C- 

t) 

1 

If  normalized,  so  that  x(r,n)  =  1  and  x(0.0)  ~  1  then 


lx(r,n)l  =  ^1  -  ^Sm^—  where  R  =  tin  (l  -  , 


which  indicates  that  zeros  of  the  ambiguity  function  will  occur  at  bins  where  n(l  -  r/T)rr  is 
an  integral  multiple  of  it. 

In  the  particular  REVGEN  case,  where  eight  updates  per  pulse  length  are  used, 

t__  mAt  _  32m  _  m 
T  T  256  8 


and  every  eighth  Doppler  bin  is  filled  (i.e.,  where  bin  n  is  a  multiple  of  8:  n  -8r),  substitu¬ 
tion  yields: 

n  ( 1  -  y)  *  =  8r  (l  =  r(8  -m)ir  , 

an  integral  value  of  tt  whose  sin  [r(8-m)ir]  =  0. 

Therefore,  in  the  special  case  of  REVGEN  where  the  Doppler  bin  separation  is  an 
integral  multiple  of  the  number  of  updates  per  pulse  length,  the  bin  centered  value  is  the 
sole  contributor  in  each  row  to  the  output  of  the  matched  filter  processor. 
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